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ABSTRACT 

Context. The Orion-KL region contains the closest examples of high-mass accretion disk candidates. Studying their properties is an 
essential step in studying high-mass star formation. 

Aims. Resolving at high spatial and spectral resolution the molecular line emission in the immediate environment of the exciting 
sources to infer the physical properties of the associated gas. 

Methods. We used the CRIRES high-resolution spectrograph mounted on the VLT to study the ro-vibrational '^CO/'^CO, the PfundyS 
and H2 emission between 4.59 and 4.72 fim wavelengths toward the BN object, the disk candidate source n, and a proposed dust 
density enhancement IRC3. 

Results. We detected CO absorption and emission features toward all three targets. Toward the BN object, the data partly confirm 
the results obtained more than 25 years ago by Scoville et al., however, we also identify several new features. While the blue-shifted 
absorption is likely due to outflowing gas, toward the BN object we detect CO in emission extending in diameter to ~ 3300 AU 
with a velocity structure close to the vis, . Although at the observational spectral resolution limit, the '^CO line width of that feature 
increases with energy levels, consistent with a disk origin. If one attributes the extended CO emission also to a disk origin, its extent 
is consistent with other massive disk candidates in the literature. For source n, we also find the blue-shifted CO absorption likely 
from an outflow. However, it also exhibits a narrower range of redshifted CO absorption and adjacent weak CO emission, consistent 
with infalling motions. We do not spatially resolve the emission for source n. For both sources we conduct a Boltzmann analysis 
of the ''CO absorption features and find temperatures between 100 and 160 K, and H2 column densities of the order a few times 
10^' cm"-. The observational signatures from IRC3 are very different with only weak absorption against a much weaker continuum 
source. However, the CO emission is extended and shows wedge-like position velocity signatures consistent with jet-entrainment of 
molecular gas, potentially associated with the Orion-KL outflow system. We also present and discuss the PfundyS and Ht emission in 
the region. 

Conclusions. This analysis toward the closest high-mass disk candidates outlines the power of high spectral and spatial resolution 
mid-infrared spectroscopy to study the gas properties close to young massive stars. We will extend qualitatively similar studies to 
larger samples of high-mass young stellar objects to constrain the physical properties of the dense innermost gas structures in more 
detail also in a statistical sense. 

Key words. Stars: formation - Stars: early-type - Accretion, accretion disks - Techniques: spectroscopic - ISM: jets and outflows - 
Stars: individual: Orion-BN, Orion source n 

1. Introduction or whether they are maybe self-gravitating non-Keplerian enti- 
ties. 

Understanding the physical structure of massive accretion disks ,,„•,,• ,s , , , , 

is one of the main unsolved problems in high-mass star for- ^hi e studies at (sub)mm wavelengths are a powerful too 

■ u • J- * »u -I- e ^ e to mamly study the cold gas and dust components on spatial 

mation. Although indirect, the main line of arguments for ac- , / , , , . . t / ^ — n u„„_k 

■ cretion disks stems from massive molecular outflow obser- ^'f^^ 1000 AU e.g., C esarom et alj | 2007[) , such 

vations that identify colhmated and energetic outflows from observations are not that well suited to investigate the inner 

high-mass young s tellar objects (YSO s, e. g., Een ning et al. ^"^^ ^™ components of massive rotating structures. In con- 

l2ro0l : iBeuther et al.ll2002bl: IZhi^r^^l2005l : fc^t al.ll2007l) . ^"^^.^ ^° *f ' mid-infrai-ed specti-al lines, e.g., ro-vibrationally 

Colhmated jet-hke outflow structures are usually attributed excited CO emission lines can trace these warm gas compo- 

to massive accretion disks and magneto-centrifugal acceler- nents. However, the spectral and/or spatial resolution was mostly 

ation. Recent 2D and 3D magneto-hydrodynamical Simula- lacking because absorption features were dommating and hence 

tions of massive collapsing gas c ores also result in th e for- Prohibiting the detection of the accretion disks in emission 

mation of mass ive accretion disks dYorke & Sonnhalteii I200I Several recent studies have further demonstrated the power of 

lKrumholzetal.Ll2007L l2009). However, it is still unclear whether high-spectral and high-spatial resolution CO mid-infrai-ed spec- 
such massive disks are similar to their low-mass counterparts, Y '^ ^.^^^^ ^ 'T'^nn^ 

hence dominated bv the central Y.SO and in Kenlerian rotation Herblg Ae Stars (e.g., | GotO et al. || 200^ . Pontoppidan et al. ,| 2008| ; 
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hence dominated by the central YSO and in Keplerian rotation, 1 firferr" 

^ ^ Ivan der Plas et al.ll2009h . 



* Based on observations of the ESO program 380.C-0380(A). 



To achieve the highest angular and spectral resolution possi- 
ble, we observed some of the closest massive disk candidates in 
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Orion at a distance of 414pc jMenten et all l2007l) - Orion-BN 
(the Becklin-Neugebauer Object), source n and IRC3 - in the 
CO V = 1 -0 transitions around 4.65 yum with the CRy o^enic hi^h 
resol ution InfraRed Echelle Spectrograph (CRIRES. lKaufl et alj 
12004 ) at the VLT. Both objects are well detected at mid-infrared 
wavelengths exhibiting various kinds of disk-signatures. 

The Becklin-Neugebauer Object: Since its detection in 
the 1960s, the BN obj ect is on e of the archetypical hi g h-mas s 
YSOsJ^ecklin & Ne ugebaued. 1 19671 : iHenning et all Il990h . 
IScoville et al.l (119831) observed the source in several frequency 
settings between 2 and 5 yum and detected molecular emission 
from several CO isotopologues (fundamental and overtone emis- 
sion), and they inferred that BN exhibits an outflow/wind as 
well as a highly confined region of molecular gas at high den- 
sities and temperatures of ~3500 K. The estimated luminosity of 
the BN object is 1 - 2 x 10^ Lp co rresponding to a B0.5 main 
sequen ce star jScoville et al.L Il983l) . More recently, jjiang et al.l 
(I2005h observed BN in polarized near-infrared emission, and 
they also identified signatures caused by an embedded accre- 
tion disk. The BN object has a high velocity along the line of 
sight of ~21 km s~^ cornpared with the cloud velocity of around 
5kms"' (IScoville et al.L[T993h . This is consiste nt with the mea- 
sured high proper motions of that object (e.g.. iPlambeck et"aLl 
Il995h . Whether the BN object is expelled from the Trapezium 
system or during a disintegration of a bound system once con- 
taining source I, so urce n and the BN object itself is still a mat- 
ter of debate (e.g., 11^120041 : lOomez et all 120051: IZapata et all 
l2009h . 

Source n: Bas ed on a bipolar radio m orphology and H2O 
maser association, iMenten & Re"idl (Il995l) suggested that this 
source may be one of the driving sources of the powerful molec- 
ular outflows within Orion-KL. Extended mid-infrared emission 
was observed perpendicular to th e outflow a xis (jGreenhill et'al], 
120041 : IShuping et all |2004|) . and iLuhmanI (|2000|) detected CO 
overtone emission. Both features are interpreted as likely be- 
ing due to an accretion disk. Source n is believed to be in an 
evolutionary younger stage than the BN object, and the lumi- 
nosity is estimated to be lower as well, of the order 2000 
(iGreenhill et al.Ll200l . 

IRC3: The source n observations s erendipitously covered 
the extended infrared source IRC3 (e.g.. lDougados etal.lll993h 
which we present here as well. At 3.6 wavelen gths, IRC3 is 
elongated in the northeast-southwest direction (Dou gados et al.L 
1199 3) and shows highly polarized near- to mid-infrared emis- 
sion (Minchin et al., 1991). The observations are consistent 
with IRC3 being a dust density enhancem ent reprocessing light 
from another source, potentially IRC2 (jPownes et al.L Il98ll : 
iMinchin et aLLll99ll : lDougados et all 1 1991 ! 



Figure[T]gives an overview of the region marking the sources 
discussed in the paper as well as the slit orientations (see also 
section|2]i. The nominal absolute positions for the three sources 
are listed in Table [1] 



Table 1. Source positions (from lPougados et al.lll993|) 



Source 


R.A. 


Dec. 




(J2000.0) 


(J2000.0) 


BN object 


05h35ml4.12s 


-05d22m22.9s 


Source n 


05h35ml4.35s 


-05d22m32.9s 


IRC3 


05h35ml 3.90s 


-05d22m30.0s 




-30 



-20 



20 



30 



-10 10 

ARC SECONDS 

REFERENCE POS. R.A. 05''35'"l4'l2 DEC -05°22'22'.'9 (J2000) 

Fig. 1. Overview image of the region . The grey-scale /T, band 
image is taken from UKIDSS DR5 (Lawr ence etall l2007l) . the 
white stellar peak positions are due to saturation. The slits and 
important sources are marked. The cross denotes the position of 
IRc2. 



2. Observations 

We obtained high-resolution spectra between 4.6 and 4.7 pm 
with CRIRES (Kiiufl et al., 2004)) mounted on UTl at the VLT 
on Paranal, Chile. Two grating settings were selected (12/-l/n, 
/lret=4662.1 and 12/-l/i, /li-ef =4676.1) to observe the spectral in- 
terval covering the '^CO v = 1 - [P(l)-P(5)/R(0)-R(8)] lines 
without gaps. 

One slit covered the BN object with a position angle of 126 
degrees east of north, whereas the second slit included source 
n and IRC3 with a position angle of 1 10 degrees east of north. 
For the strong source BN, only 40 sees (DIT=2 sees, NDIT=10) 
were required. For the second object with the weaker sources we 
had a total on-slit integration time of 20 minutes (DIT=10 sees, 
NDIT=2). For the BN observation a slit width of 0.2" was used 
while for the observation of source n a slit with of 0.4" was se- 
lected which correspond to spectral resolving powers /l/A/l of 
100,000 and 50,000, respectively. The non-AO mode was ap- 
plied since no natural guide star is available in the close environ- 
ment. The infrared seeing measured from the spectra was 0.35" 
during the BN observation and 0.45" for the source n observa- 
tion. The nod-throw of all the observations was set to 10". To 
correct for the telluric absorption lines, attached to every science 
observation a telluric standard star (HR 1666 with spectral type 
A3III) was observed. 

As extended emission in the CO lines was present in the ob- 
servations, we corrected the frames for distortion in order to get 
a wavelength solution valid for the whole chip. Firstly, the chips 
are slightly rotated with respect to the slit (ranging from 0.05 de- 
gree for chip 3 to 0.45 degree in the case of chip 4). We measured 
the position of the brightest object as function of wavelength and 
calculate from the displacement in position the rotation angle. 
Secondly, after the rotation angle has been corrected, the curva- 
ture of the slit is corrected by measuring the position (central 
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wavelength) of a sky-line as function of the spatial coordinate. 
This could be described with a 2nd degree polynomial. Using 
the IDL routines polywarp and poly_2d, the distortion was cor- 
rected. 

After the distortion correction, the raw files are processed 
by the ESO CRIRES pipeline (version 1.10.1) in combination 
with the Gasgano software. The data are dark subtracted, flat 
field corrected as well as corrected for non-linearity. The wave- 
length calibration is done using the telluric e mission lines in 
comb ination with a HITRAN model spectrum (iRothman et al.L 
|2005|) . A cross-correlation of the spectra with the HITRAN 
spectra showed that the wavelength accuracy of the spectra is 
+0.5kms-i. 

Two absorption lines are present in the standard star: the 
Pfund j6 and the Humphreys e line. Due to the strong telluric 
absorption, these lines are not trivial to remove. Therefore, we 
first reduced the spectrum without correcting for the intrinsic 
absorption lines of the standard star In the final reduced spec- 
trum these lines become eminent as emission Unes free from 
contamination by the atmosphere. The Pfund /3 line is also seen 
in our science object. However, the Humphreys e line is not 
present in the science spectrum of BN before division by the 
standard star. Therefore, it can be used to correct for the ab- 
sorption lines of the standard star We used a high resolution 
Kurucz model spectrum from an AOV star (http://kurucz.har 
vard.edu/stars.html) and scaled and shifted the model spectrum 
such that the Humphreys e profile would fit that of the observed 
Humphreys e profile in the standard star Assuming that the 
Pfund/? line scales in the same way as the Humphreys e line 
does, we divided with the model spectrum to remove the line 
contamination. 

The spectra were corrected for the earth velocity to the local 
standard of rest using rvcorrect in IRAF. The velocity correc- 
tions applied for the two observing dates were -3.8 km s ' for the 
BN data (observed on 21st October 2007) and +43.2 km s"' for 
the source n/IRC3 observations (taken on 21st February 2008). 
The velocity relative to the local stan dard of rest Vkr of O rion 
varies between 2.5 and 9kms"' (e.g.. IComito et al.ll2005h . and 
we adopt the approximate value of +5 km s~'. 

3. Results 

Toward all three sources we detected the whole suite of '^CO v - 
1 -0 lines present in the spectral window, the '^CO v = 1 -0 lines 
from R(6) to R(13) that were not blended by the '^CO v = 1 - 
lines, as well as the Pfund/? line. Figure |2]presents the complete 
spectrum toward the BN object, and Table|2]gives an overview of 
the covered lines, their wavelengths A and the lower-level energy 
state of the transitions (Ziiower/^)- In total, this setup covers a 
broad range of energy levels extending up to 504 K. 

3.1. The BN object 

3.1.1. CO absorption and emission 

Figure [3] presents a zoom compilation of '^CO data from the 
R(l) to R(8) line covering lower energy levels between 5.5 and 
199 K. If one ignores the telluric line feature at ~3kms ' two 
broad absorption features can be identified at approximately -14 
and -1-8 km s"'. Since all '^CO lines are saturated, their peak 
absorption velocities are unreliable, and we refer to the '^'CO 
data (Fig. |4]i. The peak velocity of the blue-shifted component 
is ~ -15kms"' extending from ~ -28 to ~ -8kms"'. The 
second absorption feature has its peak at ~ 8kms"', close to 



Table 2. Observed lines 



Line 


A 


^lower/^ 




jim 


K 


'^CO 


R(9) 


4.5876 


249 


R(8) 


4.5950 


199 


R(7) 


4.6024 


155 


R(6) 


4.6090 


116 


R(5) 


4.6176 


83 


R(4) 


4.6254 


55 


R(3) 


4.6333 


33 


R(2) 


4.6412 


17 


R(l) 


4.6493 


5.5 


R(0) 


4.6575 





P(l) 


4.6742 


5.5 


P(2) 


4.6826 


17 


P(3) 


4.6912 


33 


P(4) 


4.7000 


55 


P(5) 


4.7088 


83 


"CO 


R(13) 


4.6641 


504 


R(12) 


4.6741 


432 


R(ll) 


4.6782 


365 


R(10) 


4.6853 


304 


R(9) 


4.6926 


249 


R(8) 


4.7000 


199 


R(7) 


4.7075 


155 


R(6) 


4.7150 


116 


Pfund yS 


4.6538 


157954° 


H, 0-0S9 


4.6947 


7198° 



The temperatures listed for the Pfund/? and H2 - 0S9 lines do not 
indicate gas temperatures since these lines are not excited by collisions 
within thermal gas, but rather by ionizing photons and within shocked 
regions, respectively. 



BN-object. '^CO R{8) to R(l) top to bottom 




Velocity [km/s] 



Fig. 3. The '^CO R(8) to R(l) lines from top to bottom toward 
the BN object. Telluric and visi velocities are marked. The feature 
between and 8 km s"' is a telluric artifact. 



the visi of the cloud, and extends from ~ -8 to ~ 18kms"^ 
Furthermore, red-shifted from the absorption we clearly iden- 
tify a CO emission peaking in '^CO and '^CO at ~ 20kms ' 
and extending from ~15 to ~30km s The overall extent of the 
'^CO absorption and emission is from ~ -30 to — [-30kms"'. 
It should be noted that while the of the dilferent cloud com- 
ponents for the Orion-KL region var y between approximately 
3 and 9kms IScoville et all d 19831) inferred that the corre- 
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Fig. 2. CRIRES observations of the R and P CO line series around 4.65 /im. The strong absorption features are the CO lines, and 
the emission line at ~4.654/im is PfundyS. 



sponding velocity of the BN object is significantly larger around 
21k ms~- ' (con s istent with the diffe rent BN ejection scenarios, 
e.g.. lT^l2004l : iGomez et al.ll2005h . While the absorption fea- 
tures stem from the warmer protostellar envelopes and the sur- 
rounding cloud with usual temperatures of the order lOOK (see 
also Boltzmann analysis below), the ro-vibrational lines in emis- 
sion can be caused by different processes. For example, fluo- 
rescence via UV photons or resonance scattering from strong 
infrared fields c an excite these lin es without si gnificantly heat - 
ing the gas (e.g. jBlake & Boogertll20b4 : Rvde & Schoiei200lh . 
Alternatively, the ro-vibrational lines could be caused by hot- 
ter gas components (see discussion in section We note that 
the critical densities of these lines are of order lO'^^ cm"^ which 
practically implies that extended gas components can be hardly 
responsible for the emission. 

Although the blue-shifted part of the spectrum with respect 
to the Visr of the molecular cloud is slightly broader than the red- 
shifted part, nevertheless we clearly identify red-shifted absorp- 
tion as well. To first order, the blue-shifted gas seen in absorption 
can be identified with outflowing gas from the region, whereas 
the red-shifted features belong to gas infalling in the direction of 
the central source. The outflowing gas with a maximum velocity 
relative to the Vk, of ~ 35 km s"' is consistent with outflow wings 
often observed at mm wav elengths from young massive star- 
forming regions (e.g., Beu ther et alji2002bf) . It should be noted 
that the even broader outflow wings observed at mm wavelength 



toward Orion-KL exceeding +50kms"^ (e.g.. lChemin & Wrightl 
1996) are likely not related to the BN object but rather to one 
or more sources ab out 10" south-east of BN (source I, source n 
and/or SMAl, e.g., Greenhil l et al.l2004l:I.Tiang etani2005l:lBanvl 
2008 ; Beuther & Nissen 2008|)! 

How do these gene r al fea tures compare with the data pub- 
lished bv lScoville et al.l (Il983h which were observed during sev- 
eral observing runs between December 1977 and February 1981, 
hence about 30 years prior to our observations. The general CO 
line structure with two strong absorption features plus one red- 
shifted emission peak are largely the same. Also the overall ex- 
tent of the emission is quite similar Compared to our measured 
values of ~ -15, ~ +9 and ~ 20kms ' for the three compo- 
nents, IScoville et alj ([1983') report for the corresponding fea- 
tures velocities of -18, +9 and -1-20 km s"'. While two veloci- 
ties agree well, the most blue-shifted absorption peak appears to 
have shifted a little bit between the two observations. However, 
given that their spectral resolution was more than a factor 2 lower 
than that of the new CRIRES data (7 versus 3 km s ~'), we refrain 
from further interpretation of this. Furthermore, IScoviUe et al.l 
1 198 31) identify two more absorption features, one at -3kms"' 
and one at -i-30kms Regarding the -3kms"^ component we 
are not able to infer any changes because that features lies very 
close to the telluric emission which obscures any reliable sig- 
nature there. However, Fig. |3] shows that we do not detect any 
additional absorption feature blue-shifted from the 26.6 km s ' 




emission. Therefore , the +30kms ' absorption dip reported by 
IScoville et alJ (Il983h was either a transient feature or not signif- 
icant with respect to the signal-to-noise ratio. 

Since the '^CO data are so strongly saturated, for the fol- 
lowing analysis we use the corresponding '^CO data covering 
the R(6) to R(13) transitions (Fig.|2]l. Since the absorption depth 
///() is related to the optical depth r via ^ = e"'" we can directly 

estimate the optical depth of the '^CO absorption lines shown in 
Fig. m if the lines are spectrally resolved. Since the full width at 
zero intensity (FWZI) is of order 20kms"' (see above), this cri- 
terion is fulfilled with our spectral resolution AjhA - 10^ which 
corresponds to a velocity resolution of 3kms '. Except of the 
lowest '^CO R(6) line component at +9 km s ' all other observed 
'^CO absorption features have optical depths below 1. This al- 
lows us to estimate rotational temperatures via Boltzmann plots 
fr om the eq uivalent line width following the approach outlined 
in IScoville et alj (Il983l) also adopting their finite optical depth 
corrections. Using their equation (A8), the column density Ni of 
the lower-level energy state is 



Table 3. '^CO values for Boltzmann plots 



FA 

Ni = 1.13 X 10'^ — cm 

flu 



-2 



with 



F, A and /„; deno ting the op tical depth correction fol- 



s op t: 

lowing the appendix in Scoville et al.l Gl983). the equivalent line 
width (in units of cm ') and absorption oscillator strength, re- 
spectively. Because the absorption feature at ~9kms"' is af- 
fected by the telluric line, in the following we only work with 
the data from the ~ -15 km s"' line. The parameters and the cal- 
culated column densities are listed in Table [3] In a Boltzmann 
distribution, the total column density A^tot('^CO) is related to the 
temperature and the lower state column density via: 



In 



In 



jVtot(''CO) 
Q(T) 



1 £io 



with gi the statistical weight of the lower-level transition and 
Q{T) the partition function at the given temperature. 

Figure |5] presents the corresponding Boltzmann plot where 
the lower-state column density Ni divided by the statistical 
weight is plotted against the lower-state energy iiiowei divided 
by k. A linear fit to the data gives the rotational temperature Trot 





^lower/^ 


A 




F 


8.1 


Ni("CO) 




K 


km/s 


xlO-^ 






xlO'^cm-2 


BN@-15kms-' 


R(6) 


116 


3.70 


6.0 


1.18 


13 


6.0 


R(7) 


155 


3.16 


5.9 


1.14 


15 


5.0 


R(9) 


249 


1.90 


5.9 


1.08 


19 


2.9 


R(10) 


304 


1.36 


5.8 


1.06 


21 


2.0 


R(ll) 


365 


0.69 


5.8 


1.04 


23 


1.0 


R(12) 


432 


0.51 


5.8 


1.02 


25 


0.7 


R(13) 


504 


0.27 


5.8 


1.02 


27 


0.4 


source n 


@-7kms"' 


comp. 










R(7) 


155 


4.1 


5.9 


1.14 


15 


6.5 


R(9) 


249 


2.7 


5.9 


1.08 


19 


4.1 


R(10) 


304 


2.0 


5.8 


1.08 


21 


3.1 


R(12) 


432 


1.3 


5.8 


1.06 


25 


2.0 


R(13) 


504 


1.5 


5.8 


1.05 


27 


2.2 


source n 


@+5kms" 


' comp. 










R(7) 


155 


4.6 


5.9 


1.23 


15 


7.8 


R(9) 


249 


2.7 


5.9 


1.16 


19 


4.4 


R(10) 


304 


1.8 


5.8 


1.12 


21 


2.8 


R(12) 


432 


0.6 


5.8 


1.04 


25 


0.9 


R(13) 


504 


0.1 


5.8 


1.00 


27 


0.1 



The Table entries are explained in the main text. 

and the column density of '^'CO divided by the partition func- 
tion QiT). For the -15kms~' '^CO component the fitted Tjot is 
-112 + 20 K. IScoville et al.1 (1198 3^ calculated the rotation tem- 
perature for the more blue-shifted absorption component from 
the '^CO V = 2-0 transitions, and they find ~150K there, which 
is approximately consistent with our new determination. The 
other fitted pai-ameter is Mot('^CO)/Q(T) ^ 1.35 x lO'^cm^^. 
For 2-atomic molecules like '^CO the partition function can be 
approximated by Q(T = 112K) ~ kT/(hB) ~ 42 (where B is 
the rotation constant), and we get a '^CO column density of 
~ 5.7 X 10''' cm"^. Usin g furthermore the ' ^CO to '^CO iso- 
topologic ratio of 69 (e.g.7 lSheffer et al.ir2007l) . we derive a total 
CO column density of the -15 km s ' component toward Orion- 
BN of ~ 3.9 X 10"* cm"^. These column density estimates are 
in exc ellent agreement with the results derived bv lScoviUe et al.l 
d 19831) from the '^CO v = 2 - lines. 
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Fig. 5. Boltzmann plot for the '^CO v=l-0 lines toward the 
BN object. The x-axis shows the lower-level energies of the tran- 
sitions and the y-axis presents the natural logarithm of the cor- 
responding column densities divided by their statistical weights. 



How do these values compare to other observation s? The 
BN object was detected at 1.3 mm wavelength by Bla ke et al.l 
( 1199 6) at a 0.1 5 Jy level with a spatia l resolution of 1.5" x 
1.0". Following IPlambeck et all (Il995h . about 50mJy of that 
flux can be attributed to circumstellar dust emission. Assuming 
optically thin dust emission, an average dust temperature of 
100 K, a dust opacity index y6 of 2 and a standard gas-to-dust 
ratio of 100, we can calculate the total H2 column density 
(Hild ebrand, 1983; Ossenkopf & Henning, 1994; Henning et al., 
1 19951 : iBeuther et all l2002aL l2005al : iDraine et al.L l2007h . the de- 
rived H2 column density and mass within their 1 .5" X 1 .0" synte- 
sized beam are then ~ 1 .3 x 10^"^ cm"^ and 0.1 M0, respectively. 
For comparison, using a standard CO-to-H2 ratio of 8 x 10"^, the 
H2 column density estimated here from the CO IR spectroscopy 
data is ~ 4.9 x 10^^ cm"^. While the mm continuum data trace 
all velocity components along the line of sight, the near-infrared 
data trace selected velocity components. Judging from Figures 
[3] and m where we see more than 1 velocity component, the to- 
tal column density traced by the CO data is more than a factor 2 
higher. Furthermore, the CO absorption only traces the gas in the 
foreground of the near-infrared source reducing the traced gas by 
another factor 2. While both approaches are affected by system- 
atics - e.g., the mm derived column densities can be wrong by 
a factor 5 depending on the assumptions on the dust properties 
and temperatures - other reasons are more important for some 
of the differences. In particular, the mm continuum emission is 
sensitive to the cold and warm dust emission whereas the near- 
infrared absorption of the '^CO transitions useable for our anal- 
ysis traces mainly the warmer gas. Therefore, these data indicate 
that a large fraction of the gas is at relatively low temperatures. 

3.1.2. Potential signatures from tine BN disk? 

While most of the observed emission is spatially unresolved, we 
find weak extended '^CO emission toward the BN object. Figure 
|6] presents a position-velocity diagram along the slit direction 
(PA of 126° east of north, see sec. B which corresponds to the 
proposed disk orientation jJiang et al.Ll2005h . We find CO emis- 
sion around the rest velocity of the BN object of ~ 21 km s"' ex- 
tending approximately +4" in both directions. Since we do not 
identify any strong velocity dispersion, this extended emission 
is unlikely to be due to an outflow, but it may potentially come 
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Fig. 6. Position velocity diagram of the I^CO (RO) line along the 
slit direction (PA of 126° east of north, see sec. |2]i plotted in 
logarithmic intensity scale. Around offset the figure is domi- 
nated by the continuum emission, but outside ~ +2" the emis- 
sion stems from the CO line. Positive offsets go in northwest 
direction. 



from the proposed disk (Ijiang et alil2005h . The measured extent 
of ~ 8" would then correspond at the given distance of Orion of 
414 pc to an approximate disk diameter of ~ 3300 AU or a disk 
radius of 1650 AU. While such disk size would be comparably 
la rge with respect to typ ical low-mass disks (e.g., several reviews 
in lReipurth et al.ll2007h . it is consistent with measured sizes of 
rotating structures in other high-mass star-forming regions (e.g. , 
Schrever et al. '2002'; 'C esaroni et all l2005l : iBeltrin et alJ l2006l : 
Beuther & Walsh 2008). 

While we do not identify extended emission is the rarer 
'^CO isotopologue, the CO & '^CO emission feature around 
20kms ' exhibits an additional interesting feature in the '^CO 
data (Fig. HJi: While the emission feature is at the edge of the 
bandpass for the '^CO R(6) line, it remains undetected for the 
R(7) line and only comes up for the lines greater R(9). This 
emission feature is also visible in all '^CO lines (R(l) to R(8)) 
shown in Figure [3] but it does not exhibit any significant shape 
change there. Likely, the '^CO emission is optically thick and 
traces only an outer envelope that does not show big variations 
with excitation temperature. In contrast to that, the measured 
Gaussian line width of that component for the more optically 
thin i^CO R(9), R(10), R(ll), R(12) and R(13) lines are 4.3, 
5.3, 6.5, 6.4 and 6.6 km s"', respectively. Although errors on the 
line width are difficult to quantify because the emission is at the 
flank of the strong absorption feature, and furthermore our spec- 
tral resolution is only 3kms"', the data are indicative of a line 
width increases with excitation temperature of the line that ap- 
pears to saturate for the highest detectable transitions R(ll) to 
R(13). While for a centrifugally supported disk with a central 
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Fig. 7. Hydrogen recombination Pfund/? line toward the BN ob- 
ject. The histogram shows the data, and the three other 
lines present a two-component Gaussian fit to the data. Both 
Gaussians are shown separately as well as the combined fit. The 
FWHM Av of both components is given in the figure. The big 
dots reproduce the best fit to the Bra line (their Figure 7 re- 
scaled to our normalized Pfund^S spectr um) for an optically t hin 
outflow with a velocity law v oc r^^^^ by lScoville et al.l (Il983h . 

mass of ~10Mo the circular velocity at 3300 AU is relatively 
small of order ~ 1.6km s ', the measured line width increase is 
consistent with a rotating disk structure because the inner region 
with larger rotation velocities should have higher temperatures 
caused by the exciting central star. One has to keep in mind that 
for purely thermal excitation, rotation would not cause such an 
effect because the upper levels - all in the 3000 K regime - do 
not exhibit big relative excitation differences. Hence, they do not 
trace significantly different regions of a disk then. In contrast to 
that, for lines that are excited by UV fluorescence such an ef- 
fect is possible because the excitation mechanism via an elec- 
tronic excited state tends to preserve the level population of the 
V = rotational states (e.g., Brittain et al. 2009). Therefore, in 
this case the R(9) to R(13) lines are sensitive to gas temperatures 
in a relatively broader range between 150 and 500 K (Table |2|i. 
An exact reproduction of the line width increase would require a 
detailed disk model, including its density and temperature struc- 
ture. Since we cannot constrain these parameters well from our 
data, this is beyond the scope of this paper. Although the line 
width increase is below our nominal spectral resolution element, 
qualitatively the observations are consistent with a rotating disk 
structure. 

3.1 .3. The PfundyS line toward the BN object and H2 
emission nearby BN 

Figure I2] presents a zoom into the PfundyS line. The line shape 
consists of two components, one central Gaussian component 
and broad line wings. It is possible to fit the whole profile rela- 
tively well with a two-component Gaussian fit where the central 
component has a FWHM Av = 39.3 km s ' and the broad com- 
ponent has a FWHM of Av - 175.5 km s ' . The full width down 
to zero intensity is approximately 340km s"' (between -170 and 
170 km s '). Close to the peak of the profile at around 25 kms"', 
the spectrum exhibits a small dip which is likely an artifact from 
the telluric corrections (see section|2ll. Considering this, the peak 
of the central Gaussian fit at ~16.5 km s ' is still consistent with 
the velocity derived for the BN object by IScoville et al.l (Il983h 
of ~21 kms"'. 

Figure [71 a l so pre sents as thick dots the best fit obtained by 
IScoville et al.l (Il983l) for the Bro- line (their Figure 7 re-scaled 



-80 -60 -40 -20 20 

Velocity (km/s) 

Fig. 8. Position velocity digram of the H2 emission ~ 4.3" south- 
east of the BN-object a ssociated with the infrared source IRC15 
(iDougados et al.Lri993h . 

to our normalized spectrum). Their model consists of a super- 
sonic, optically thin decelerating outflow with a velocity law of 
V oc r^^/^. It is remarkable how well the shape of their Bra line 
obtained ~ 30 years ago corresponds to the shape of the newly 
observed PfundyS line. Hence, these Pfund/? data are also con- 
sistent with their outflow model. 

Furthermore, approximately 4" south-east of the BN object 
we detect H2 emission from the H2 - 0S9 transition with 
Eiovia /k = 7198 K. This H2 emission feature is spatially associ- 
ated with the c orresponding H2 knot in near-infrared H2 images 
(e.g.. lNissen e t al. 2007 ) as well as wit h the mid-infrared source 
IRcl5 reported by Shuping et al.l (1200 4). Figure [8] shows a po- 
sition velocity diagram of this feature. While we do not resolve 
any spatial structure of the H2 knot, it shows a very broad veloc- 
ity extent going to blue-shifted velocities of ~ -80kms"', in ex- 
cess of the CO absorption features measured toward the BN ob- 
ject. Since we are mainly interested in the BN object, source 
n and IRC3, we refrain from further analysis of this offset H2 
emission. 

3.2. Source n 

3.2.1. CO absorption and emission 

Due to the different observing dates, for source n (and IRC3 
in the following section) the telluric lines are almost entirely 
shifted out of the CO spectrum. Figures |9] and [10] present the 
corresponding '^CO R-series lines and the detected '^CO lines 
with their associated optical depths. In comparison to the previ- 
ous BN data, for source n we only identify one broad absorption 
feature in the ^^CO data which is dominated by blue-shifted out- 
flowing gas. It is interesting to note that the most blue-shifted 
absorption does not have a Gaussian shape but rather a more 
extended wing-like structure. In the framework of accelerated 
winds with distance from the driving star or disk, such a spec- 
tral behavior would be expected. These accelerated winds in- 
crease in velocity with distance from the star Simultaneously, 
with increasing distance from the driving source the density of 
the surrounding gas and dust envelope also decreases, lowering 
the corresponding absorption depth at those velocities. Hence, 
in this picture, higher-velocity gas should show shallower ab - 
sorption features in the spectra (e.g.. Lame rs & Cassinellil 19991) . 
Furthermore, we also identify a red-shifted emission component 
at — hl5 kms However, compared to the BN object, where the 
emission feature can be identified in all transitions, for source 
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Fig. 9. The ^^CO R(8) to R(l) lines from top to bottom to- 
ward source-n. Telluric and visi velocities are marked. The linear 
slopes and dips around -40kms ' associated with the R(8), R(6) 
and R(4) lines (for the latter extending from ~-5 to ~-65 km s"') 
are artifacts due to the telluric corrections. 

n it is more prominent in the higher excited lines. The total 
width of the CO absorption and emission is ~90kms ' rang- 
ing from approximately — 65 to ~-H25kms"' (broader than for 
BN). The blue-shifted end is less well determined because of tel- 
luric line contamination. Nevertheless the blue-shifted outflow 
part of the spectrum extends ~ 70kms"' from the assumed visi- 
of -1-5 km s"'. This value exceeds that measured toward BN by 
about 20 km s"' . Although at the edge of the telluric line contam- 
ination, we tentatively identify an additional discrete absorption 
feature at — 35 km s"' . 

Since the '^CO line is again always saturated (Fig. |9]l, for 
the additional analysis we use the corresponding '^CO data 
(Fig. [Tol l. Because of line-blending and telluric emission, we 
only measure five '^CO lines for source n. The broad, saturated 
'^CO absorption feature is clearly resolved into two compo- 
nents with approximate velocities at — 7 km s"' and — h5 km s"' . 
While the ~-H5kms"' shows the deeper absorption features for 
the lower energy R-lines, it is interesting to note that the single 
absorption peak observed for the '^CO R(13) line is associated 
with the ~-7kms"' component, indicating hotter gas at more 
blue-shifted velocities. Following the approach outlined for BN 
in section U. 1.11 the derived optical depth of '^CO is ~0.6 at its 
highest. 

Since the telluric lines are far away in velocity space (see 
Fig. |9]l, for source n we can conduct the Boltzmann analysis for 
both absorption components separately. The measured equiva- 
lent widths A for both components are listed in Table I^Tl. Again 
calculating the '^^CO v=l column densities (Table |3ll and pro- 
ducing Boltzmann plots (Fig. fTTT i. we find that for both veloc- 
ity components fits to all 5 data points are less good than in 
the case of the BN object. This is likely due to the fact that fit- 
ting 2 Gaussians to the broad R(13) line, which does not show 
two well separated absorption features anymore, may overesti- 
mate the contribution of the ~-7kms ' component and hence 
underestimate the contribution from the — 1-5 km s ' component. 
Therefore, we also fit only the 4 lower-energy transitions be- 
tween R(7) and R(12) resulting in better fits (Fig. fTTT i. As ex- 

' The equivalent width A is measured by simultaneous Gaussian fits 
to both absorption features. 



pected from the diff'erent behavior of the two absorption com- 
ponents with increasing energy levels, the derived rotation tem- 
peratures using the 4 lower energy lines of the ~-7kms"' and 
~H-5kms"' are ~163+20K and ~103+10K, respectively. As a 
comparison, rotational temperatures measured at submm wave- 
length from CH^ OH emissio n lines toward source n are around 
200 K dBeuther et al.L l2005b^ . Considering that these measure- 
ments are conducted with diff'erent molecules and very different 
observational techniques, the overall range of similar tempera- 
tures derived at mid-infrared and submm wavelengths is reas- 
suring for the complementarity of such multi-wavelengths ob- 
servations. 

Fitting only the 4 lower-energy lines, the column den- 
sity of '''CO divided by the partition function QiT) results in 
Mot('^^CO)/Q(T) * 1.0 X lO'^cm-^ and Mot(13CO)/Q(T) * 
2.5 X 10'^ cm"^ for the ~-7kms"' and — 1-5 km s"' compo- 
nents, respectively. Approximating again the partition function 
by Q{T) ~ ktlQiB) (section [TTB we have e(163K) * 62 and 
Q(103K ) ^ 39, and using the '^CO to '^CO isotopologic ratio of 
69 (e.g., Sheffer et ai]|2007l') . the derived total CO column densi- 
ties toward source n are 4.3 x 10'^ cm"^ and ^ 6.7 x 10'** cm"^, 
of the same order as for the BN object. With the CO-to-Ha con- 
version factor of 8 x 10"^, the corresponding H2 column densities 
are 5.3x10^^ cm"^ and 8.4x10^^ cm"^. Adding these two column 
density values, the total H2 column density traced toward source 
n by this near-infrared observations is ~ 1.4 x lO^"* cm"^, more 
than an order of magnitude below the H2 column densities de- 
rived with the Submillimeter Array at 865 jjm (~ 6 x 10^"^ cm"^, 
Beuther et al. 2004). Similar to the BN-case discussed in sec- 
tion l3.1.n this shows that the near-infrared data mainly trace the 
warmer gas whereas the (sub)mm continuum observations trace 
warm and cold gas components. 

3.2.2. The PfundyS line 

The PfundyS line is also detected toward source n (Fig.fTSI). and 
we can fit a Gaussian to the recombination line with FWHM of 
Av X 24.3 km s"', a width down to intensity of » 70kms"' 
(from ~ -35 to ~ +35kms"') and a central velocity of 
~0 km s ' . The line is covered by 2 CRIRES chips, and while the 
general PfundyS line shape is the same for both chips, the central 
dip at ~0 km s ' cannot independently be reproduced. Therefore, 
the dip is likely only an artifact due to insufficient signal to noise. 

The thermal line width of a lO'^K hydrogen gas is 
~21.4kms"'. Convolving that with the 6kms"' spectral resolu- 
tion, the observable thermal line width should be ~22.2kms 
Therefore, the measured line width does not exceed much the 
thermal line width of an Hii region. Hence, the Pfund /? emission 
from source n does not exhibit strong signatures from a wind but 
is rather consistent with a thermal Hii region. We also cannot ex- 
clude that the Pfund yS emission toward source n is contaminated 
by more broadly distributed Orion nebula emission. 

3.3. IRC3 

As outlined in the Introduction, in contrast to the BN object and 
source n, the source IRC3 is unlikely to be a YS O, and it is 
not prominent in any typi cal hot core tracer (e.g., iBlake et "sH 
[1999 :1 Beuther et al ]|2005bh . The source rather resembles a more 
extended dust density enhancement that reprocesses light from 
other sources, potentially IRC2. Figure [T3]presents a part of the 
2D slit spectrum clearly showing the extended nature of the CO 
emission toward IRC3 in contrast to the point-like continuum 
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Fig. 10. CO optical depths toward source-n. The absorption feature at — H25kms of the R(7) line is likely an artifact due to 
telluric lines. 
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Fig. 13. Original data from the slit covering source n and IRC3. Source n is the bright continuum source at offset ~ 4" whereas 
IRC3 exhibits the extended CO emission around offset ~ 10". 



structure from source n. The continuum emission from IRC3 is 
very weak compared to that from source n. 

Figure [14] shows the CO lines extracted as an average spec- 
trum of length 5.16" toward IRC3. While we again see a broad 
absorption feature against the weaker background extending 
from ~ -5kms"' out to the telluric contamination at about 
-40km s ', IRC3 shows a broad and strong emission feature 
peaking at ~ 2kms"^ Since this emission is extended over sev- 
eral arcseconds. Figure [15] presents the position velocity dia- 
grams of selected '^CO and all '''CO lines. In contrast to the 
BN object where the extended CO emission is just around the 
rest velocity of the star (Fig. [6]), here we see a clear trend of in- 
creasing velocity with increasing distance from the continuum 
peak (the so-called Hubble law of outflows). Furthermore, in 
particular the '^CO position velocity diagrams show a twofold 
structure with one velocity increase to values > lOkms"' at off- 
sets of ~ -1.8" and another increase to similar velocities at 
~ -3.3" offset from the main continuum peak. To guide the 
eye, these wedge-like structures are sketched in the top-right 
panel of FigurefTsTl. Such multiple wedge position velocity struc- 
tures are what jet-bow-s hock entrainment models for molecu- 
lar outflows predict (e.g., lArce et aL I I2OO7I) . Furthermore, simi- 
lar to source n, the CO absorption shows also for IRC3 the ex- 
tended wing-like features toward the most blue-shifted absorp- 
tion. Similar to the pv-diagrams, where we see acceleration of 
the gas with distance from the source, these wing-like absorption 



can also be interpreted in the framework of accelera ted winds 
(see also section 13.2.11 or iLamers & C assinell? '1999). In sum- 
mary, IRC3 is not only a dust density enhancement scattering 
the light from another source, it may also be part of the fa- 
mous outflow emanating from the Orion-KL region. Based on 
these data, we cannot clearly define the driving source of the 
CO outflow feature. While the continuum peak of IRC3 could 
be driving an outflow, the CO emission may also be part of the 
large-scale outflow where different driving sources have been 
proposed for in the literature, e.g., source I, SMAl, or the dis- 
integration of a bou nd system once containing sources I, n and 
the BN-object (e.g., Menten & Reid 19951 iGomez et al]|2005l : 
Beuther & Nissen 2008; Zapata et al. 2009|). While flie position 
velocity structure of the '^CO lines largely resembles that of the 
'^CO lines (Fig. [TSl l, it is interesting to note that the strongest 
'^CO feature is at higher velocities and relatively distant from 
the continuum peak (offset ~ -2.8" and velocity ~ -)-7kms"'. 
This difference is likely caused by the varying optical depth of 
the '^CO and '^CO transitions. Since the velocity of the emission 
feature is very close to the visi of the cl oud (typically between 
5 and 9kms e.g.. IComito et al.ll2005l) . self-absorption of the 
most abundant '^CO isotopologue may veil this feature, whereas 
the less abundant '^CO can be utilised to penetrate more deeply 
into the cloud and hence detect the emission better 
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Fig. 15. Position velocity cuts along the slit axis for 1RC3 with a PA of 110 degrees east of north. The top row shows diagrams for 
selected '^CO lines, and the bottom row presents the '^CO data. For '^CO the scaling is adapted to highlight the weaker extended 
emission in contrast to the continuum at offset 0". The wedge-like outflow structures are sketched in the top-right panel. 



4. Discussion and conclusion 

High spectral resolution mid-infrared observations allow us to 
infer several important characteristics of some of the major 
sources within the Orion-KL region. The fundamental CO lines 
are largely dominated by absorption from the individual YSO 
envelopes and their surrounding gas cloud. These absorption 
features are blue- and red-shifted with respect to the vis^ of the 
molecular cloud indicating that outflowing and inflowing gas are 
simultaneously present toward our target sources. However, we 
also identify interesting emission features. Together they con- 
firm the youth of the sources where infall and likely accretion 
are still ongoing. 



BN object: For the BN object, our data confirm, at double 
the spectra l resolution, s e veral of the assessments conducted 
already by IScoville et al.l (Il983l) . However, also discrepancies 
arise. For example, an absorption feature reported previously 
around -1-30 km s"' is not found in the new data, indicating ei- 
ther transient components or poor signal-to-noise in the older 
data. From a Boltzmann analysis, the rotational temperature is 
around 112 K, and we derive CO column densities of several 
times 10'** cm"^, well in agreement with the older results based 
on '^CO v=2-0 by Scoville et al. ( 1983). Using standard CO-to- 
H2 conversion factors, these column densities are about an or- 
der of magnitude below H2 column density estimates based on 
mm continuum emission. As discussed in section 13.1.11 while 
systematics may account for some of this discrepancy, the main 
difference is that the mm continuum emission is sensitive to cold 
and warm dust, whereas the near-infrared absorption lines trace 
only the warm gas components. 



We also identify extended CO emission, likely emanating 
from the close environment of BN. While the absorption stems 
from warm gas with temperatures of the order 100 K, getting the 
ro-vibrationally lines in emission implies already much higher 
temperatures for the corresponding gas components. The veloc- 
ity range of the CO emission between 15 and 30 km s"' en- 
compasses the previously inferred veloc i ty of the BN object 
of ~21kms Following IScoville et al.l (Il983h . who covered 
more J-transitions (fundamental and overtone emission) than we 
do, the rotational temperature of the compact emitting gas is 
around 600 K, significantly higher than the above derived rota- 
tional temperature of the absorbing envelope gas. As outlined 
in section 13.1.11 high critical densities are required to produce 
this line emission by pure thermal collisional excitation, and 
other processes like UV fluorescence or resonance scattering 
from infrared emissi on may also contribute to the emission lines. 
IScoville et alj ( Il983h infer an approximate size for that emission 
of ~20 AU, whereas we now resolve the emission extending to 
~ +4", corresponding to a diameter of ~ 3300 AU. This discrep- 
ancy can be explained by the large dynamical range between the 
compact emission close to the source itself and the compara- 
bly ver};_we^more extended emission. The small size found 
bv lScovilleetan ri983) can be attributed to the strong compact 
emission which also emits in the overtone bands, whereas our 
more sensitive new data also detect the weaker extended features 
due to the high dynamic range available with CRIRES . While the 
absolute ratio of the peak emission (CO plus continuum) to the 
extended CO emission is ~210, the ratio of the continuum sub- 
tracted C O emission toward BN compared to the extended CO 
is still 70. IScoville et al.l(ll983i) were not sensitive enough to de- 
tect this faint extended emission compared to the strong central 
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Fig. 11. Boltzmann plots for the '^CO v=l-0 lines toward 
source-n for the — 7kms ' and the ~+5kms ' components in 
the upper and lower panels, respectively. The x-axis shows the 
lower-level energies of the transitions and the y-axis presents the 
natural logarithm of the corresponding column densities divided 
by their statistical weights. In both panels, the dashed line shows 
a fit to all five data points, whereas the full line presents a fit to 
only the four lower transitions. 
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Fig. 12. Hydrogen recombination Pfund/? line of source n. 



source. Although at the spectral resolution limit of our obser- 
vations, the measured line width increase of the '^CO emission 
feature with increasing excitation temperature is consistent with 
a disk origin of the CO emission. Therefore, these parameters 
make the CO emitting gas potential ly to be associate d with a 
disk around the BN object (see also iJiang et al.ll2005l) . Further 
supporting the disk-outflow interpretation for the BN object are 
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Fig. 14. The '^CO R(8) to R(l) lines (except of the R(4) Hne 
from top to bottom toward IRC3. Since the emission is extended, 
this spectrum is an average over 5. 16" along the slit. Telluric and 
visr velocities are marked. 



the Pfund /? hydrogen recombination line data which show broad 
high- velocity line wing s consisten t with the decelerating outflow 
scenario proposed bv .Scoville et al.. C1983,) . 



Source n: While bright sources like the BN o bject were al- 
ready feasible to be observed a while ago (e.g., IScoville et al.l 
1983) weaker sources like source n or IRC3 allowed reason- 
able high-spectral-resolution spectroscopy only with the advent 
of recent instruments like CRJRES on the VLT. The general pic- 
ture for source n is relatively similar A single broad absorption 
feature extending approximately to -65kms ' traces mainly 
the molecular outflow whereas red-shifted emission likely stems 
from an inner infalling and accreting envelope/disk. This pic- 
ture is consistent with disk/outflow proposals for this source de- 
duced from cm and mid-in frared wavelength imaging projects 
(e ^lMenten & Reidlll995t iGreenhiU et alj|2004l : |Shuping et alj 
I^OCMI). The clearly resolved double-peaked '"'CO structure al- 
lows to conduct the Boltzmann analysis for both components. 
We find that the colder component has higher H2 column den- 
sities (~ 103 K & ~ 8 X 10^^ cm"-) compared with the second 
warmer component (~ 163 K & ~ 5 x 10^- cm"^). As discussed 
above for the BN object, differences between the infrared CO 
derived column densities and those estimated from mm wave- 
lengths observations may arise because both tracers are sensitive 
to gas (and dust) at different temperatures. In contrast to the BN 
object, the Pfund/? emission from source n is consistent with a 
thermal Hii region without a strong wind component to the hne 
shape. 



IRC3: The observational signatures from the dust density en- 
hancement 1RC3 are very different compared to the two pre- 
viously discussed sources. The continuum from IRC3 is much 
weaker, nevertheless we detect CO absorption between ~ -5 
and ~ -40 km s"' . However, more importantly, we clearly detect 
extended CO and '''CO emission over scales of ~ 4". This ex- 
tended CO emission shows a multiple wedge-like velocity struc- 
ture consistent with jet-entrainment models of molecular out- 
flows. Hence 1RC3 may be not merely a dust density enhance- 
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ment, but it may be part of the famous outflow from the Orion- 
KL region. 

Limitations and future: One shortcoming of the data is that for 
our primary targets, the BN object and source n, we could not 
spatially resolve the inner regi on of the emissio n as done in 
the lower-mass case presented in lGoto et alj (l2006h . The reasons 
may be different for the two sources. Source n is probably still 
too young and too deeply embedded so that the envelope over- 
whelms any emission from the embedded disk itself. The BN 
object has the advantage that its velocity of rest is offset from 
that of the cloud by more than 10 km and hence absorption 
could be less of a problem for such kind of source. However, BN 
is likely significantly more evolved, and the continuum-to-line 
ratio is so high that we cannot reasonably filter out the contin- 
uum emission. Hence, we cannot well study the inner region of 
the proposed disk. Nevertheless, observing higher J-transitions 
as well as the CO o vertone emiss ion with todays higher sensitiv- 
ity compared to the lScoville et al., (1983) observations will likely 
constrain the proposed disk structure in more detail. 

How to proceed now if one wants to do similar-type fun- 
damental CO line studies of disks in high-mass star formation? 
On the one hand, it is important to not select too young sources 
because their envelopes will likely almost always "destroy" 
the emission signatures. There may exist exceptions where one 
views straight through the outflow cavity face-on toward the 
disk. On the other hand, for more evolved regions, the continuum 
emission can be very strong or maybe the remaining disk size 
can be reduced again making the spatial resolution a problem. 
Therefore, in addition to very careful target selections, adding 
AO to achieve the best spatial resolution will be a crucial element 
for such kind of studies in the coming years. Furthermore, in par- 
ticular for the mostly saturated '^CO lines, it will be important to 
extend the spectral coverage to also observe higher excited CO 
lines, that will likely not saturate anymore, as well as CO over- 
tone emission. This will allow us to better assess the hotter gas 
components and hence to conduct a more detailed analysis of 
the '^CO data themselves. On longer time-scales, the ELT with 
its proposed mid-infrared instrument METIS promises orders of 
magnitude progress in this field based on its superior sensitiv- 
ity and spatial resolution. With this instrument, we will be truly 
capable to resolve the gas signatures of accretion disks around 
(high-mass) YSOs. 
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